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Abstract 

A symmetry breaking mechanism is shown to occur in an array composed of sym- 
metric bistable Lorenz units coupled through a nearest neighbour scheme. When 
the coupling is increased, we observe the route: standing oscillating travelling 
front. In some circumstances, this route can be described in terms of a gluing of two 
cycles on the plane. In this case, the asymptotic behaviour of the velocity of the 
front is found straightforward. However, it may also happen that the gluing bifur- 
cation involves a saddle-focus fixed point. If so, front dynamics may become quite 
complex displaying several oscillating and propagating regimes, including chaotic 
(Shil'nikov-type) front propagation. Finally, we compare our results with different 
couplings and with other discrete bistable media. 
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1 Introduction 



Several areas of science use model equations of the reaction-diffusion type. Moreover, besides 
its special interest in some fields (chemistry [1,2], biology [3], ...) reaction-diffusion equations 
are considered as abstract models for pattern formation [4]. 

Real systems are frequently composed of discrete elements, material models and biological cells 
are just two examples. Therefore, it is natural to deal with the discrete space version of the 
react ion- diffusion equation. In this point, it must be emphasized that discreteness may manifest 
itself in a strong way, such that some phenomena are not simple analogies of the continuous 
ones. 
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Here, wc focus on a discrete reaction-diffusion model such that its local dynamics (reaction 
term) is bistable, with two stable fixed points. Particularly important examples of bistable 
dynamics are found in optics [5,6], chemical systems [7], and biology [3,8,9]. The main interest 
of these systems hes in the behaviour of the fronts that constitute the boundaries of the domains 
of both (stable) states. Intuitively, the front will move from the most to the less stable state, 
enlarging a domain and shrinking the other. In fact, great attention has been devoted to the 
phenomenon of 'propagation failure' or 'pinning' [10-15], because it is usual that in discrete 
systems some threshold must be surpassed to achieve propagation. 

But, although one could expect that propagation does not succeed in the case of symmetric 
bistability, it was demonstrated that such possibility exists. Hagberg and Meron [16] studied a 
continuous bistable reaction-diffusion system, the FitzHugh-Nagumo model, showing that, in 
the symmetric case, propagation occurred after a symmetry breaking front bifurcation. In this 
scenario, the front breaks its symmetry through a pitchfork bifurcation, at the same time that 
propagation is initiated. Also, an analogous bifurcation was found for the complex Ginzburg- 
Landau equation [17] with the name of nonequilibrium Ising-Bloch bifurcation. 

In a recent work, we have reported [18] the existence of front propagation in a discrete symmet- 
ric bistable medium. The transition is exclusive for a discrete system and it is possible thanks 
to the multi- variable nature of the local dynamics. One variable systems are not able to break 
symmetry [10] even if the function that describes the local dynamics is not continuous [19]. It 
was shown that the mechanism leading to propagation is not a pitchfork bifurcation. Contrast- 
ingly, it consists on a Hopf bifurcation followed by a global bifurcation, that is equivalent to 
a gluing bifurcation of cycles. In correspondence, the velocity of the front shows a logarithmic 
dependence close to the onset. The aim of the present contribution is two-fold. On one hand, 
a more detailed analysis of the transition is undertaken, explaining which is the role of the 
stationary states of the system. On the other hand, the new features of the transition when a 
saddle-focus fixed point mediates the gluing bifurcation are investigated. 

The paper is organized as follows: Section 2 presents the model. In Sec. 3 an overview of 
the transition is presented. In Sec. 4 a reduced cylindrical coordinate system is defined, and 
the mechanism of gluing of cycles underlying the transition between oscillating and travelling 
fronts is explained. In Sec. 5, several non-trivial front dynamics are shown, including the chaotic 
motion of the front. Finally, Sec. 6 is devoted to further considerations, while in Sec. 7 we present 
the conclusions. 



2 The system 

We consider a discrete reaction-diffusion equation in one dimension: 

= /(r,) + yr(r,.+i + r,_i - 2rj), (1) 
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where D accounts for the couphng strength between neighbours, and the couphng matrix F says 
which variables get coupled. In what concerns the local bistable dynamics, we deal with the well- 
known Lorenz oscillator [20] as unit cell. It exhibits several characteristic behaviours depending 
on its internal parameters [21]: monostability, bistability, limit cycle, 'butterfly chaos', 'noisy 
periodicity', etc. We consider the range where the system contains two stable symmetry related 
fixed points (bistability). The coupling matrix F is taken with all elements zero except one: 
^ — Iki — 5k2^i2- Therefore, the equations describing the temporal evolution of the array are: 

Vj = rxj - XjZj - Vj + — (yj+i + - 2yj) j^l,...,N (2) 

The standard values are chosen for the a and b parameters: a = 10 and b = 8/3. Then, the 
Lorenz system contains two stable foci C± = {±\Jb{r — 1), ±^Jb{r — 1), r — 1) for r in the range 
(1.35, th); a,t r — th ^ 24.74, C± become unstable through a subcritical Hopf bifurcation. Also, 
it is important to note that there exists a saddle fixed point located at the origin for r > 1. 

A fourth-order Runge-Kutta method, with time step 0.01, was used to integrate Eq. (2). A step- 
like initial condition is imposed to the system. Results on front propagation do not depend on 
the boundary conditions provided that the array is large enough. Thus, if one wishes the system 
evolve for long time in the regime of travelling front, one should move the boundary whilst the 
front propagates. Also, one may imposed to the system periodic boundary conditions, with two 
fronts diametrically opposed as initial conditions. If the initial condition satisfies Xj = —Xj-^.N/2, 
Uj = —yj+N/2, Zj = Zj+N/2, both fronts propagate always in the same direction, avoiding its 
mutual annihilation. 



3 Standing, oscillating and travelling fronts 

Three main front dynamics are found when varying the parameter r (r < th), and the coupling 
strength D {D > 0). Wc distinguish among standing (static), oscillating and travelling fronts. 
In Fig. 1, a centered step-like initial condition is imposed to the system with free ends and 
r = 14; two different non-static regimes are achieved depending on the value of the coupling 
strength D. Note that the Lorenz model is symmetric and the coupling that appears in Eq. (2) 
destroys this local symmetry but preserves the global reflection (or Z2) symmetry: 

{xi,yi, zi, xn, Vn, zn) — > {-xi, -yi, zi, -xn, -yN, zn) (3) 

Therefore, for large enough, both senses of propagation are equally likely. A diagram of the 
different regions on the r — D plane is shown in Fig. 2. It may be seen that static fronts are 
found for D and/or r small, whereas travelling fronts appear for large D and r. Also, it is 
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Fig. 1. Spatio-temporal evolution of the front for two different values of the coupling strength D and 
r = 14. A step-like initial condition is imposed for an open array of = 50 oscillators. When the 
coupling surpasses the critical value D = Dth ~ 39.63 the front propagates through the array shifting 
the oscillators from C+ to C_. The time interval shown is 100 time units. 

significative that oscillating fronts always exist between the standing and the travelling fronts 
regions. The line {Dos), that defines the boundary between the standing and the oscillating 
regions, approach D = Q as r ^ Th, since the nature of the fixed points {C±) deeply influences 
the dynamics of the front. This is not very surprising, but indicates that oscillating fronts will 
be usually found in bistable systems with stable foci, rather than nodes. Somehow, the stability 
properties of the fixed points are transmitted to the front. Also, it is significative that Dth (the 
threshold for traveUing fronts) seems to diverge asymptotically at r = Too ~ 13.5. 

Instead of visualizing the system as in Fig. 1, we show now which is the behaviour of the 
different units when the couphng is increased from zero. It is clear that for D = the step-like 
solution that we were considering consists of = C+, j — 1, ... ,25 and = C_, j = 26, . . . , 50 
(N — 50). As the coupling increases this solution can be smoothly continued. But there exist 
two constrains for the stationary solutions, that come from Eq. (2): 

^j=yj, ^ Vje{l,2,...,Ar} (4) 

Therefore, all the units lie on a parabola that passes through C± and zero. However, when Dos 
is reached the front undergoes a Hopf bifurcation, and all the units start to oscillate leaving the 
mentioned parabola. The projection of the parabola onto the plane {x, y) is a straight line, the 
bisectrix of the first and third quadrants. However, to better visualize oscillations out of the 
parabola, in Fig. 3 we performed a 45° rotation of the reference system. Obviously those units 
closer to the front oscillate with larger amplitude that those located far from the front that 
are almost insensitive to the bifurcation. When Dth is reached a multiple collision occurs; the 
orbits of neighbouring units collide creating two "channels" going from C± to Cip. Like in the 
symmetric FitzHugh-Nagumo model studied by Hagberg and Meron [16], the travelling front is 
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not symmetric. However, in our case the mechanism leading to propagation is not a pitchfork 
bifurcation, instead, a Hopf bifurcation that creates the oscillating front is the precursor of the 
travelling front. 

The description provided by Fig. 3 is illustrative, but it is incomplete if one does not realize 
that besides the static solution that loses its stability through a Hopf bifurcation, there must 
exist another (unstable) stationary state that mediates the multiple collision of cycles. There- 
fore, one must search the stationary solutions monotonic in {xj} and {yj}. It is not difficult 
to find out that only two monotonic solutions, called stable and unstable dislocations, exist 
(discarding spatial translations). They are the continuation of the solutions in the uncoupled 
limit: (. . . , C+, C+, C_, C_, . . .) and (. . . , C+, C+, 0, C_, C_, . . .). We shall refer to them as A- 
and -B-state, respectively, and Fig. 4 shows a sketch of them. It is the 5-state which mediates 
the transition at D^h, and it will be shown below that the stability properties of this solution 
will determine important features of the onset of the wavefronts. 
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Fig. 3. Projections onto the reference system ((xj + yi) / \/2, {—Xi + yi) / \/2) of all the oscillators of the 
array with step-like configuration and r = 14. From top to bottom, standing, oscillating, and travelling 
cases are shown. In the travelling case, the line is the trajectory followed by each oscillator going from 
C_i_ to C- as the front advances. 




Fig. 4. A sketch of the two monotonically decreasing stationary solutions for > 0. a) ^-state, 
continuation of the state (. . . , C+, C+, C_, C_, . . .) at D = 0. b) i?-state, continuation of the state 
(...,C7+,C7+,0,C_,C_,...) at Z) = 0. 
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4 Gluing of cycles 



In this section we demonstrate how the transitions presented in the previous section may be 
described in the context of a gluing bifurcation, in which two limit cycles become a two-lobed 
cycle by involving an intermediate saddle point. The gluing bifurcations is usual in systems with 
Z2 symmetry [22,23], because provided the existence of the symmetry it becomes a codimension- 
one bifurcation. 



4-1 Cylindrical coordinates 



Figure 3 shows some kind of collective motion when the front oscillates. Then, the whole set 
of variables {xj,yj,Zj) is discarded, and instead, we chose to build a reduced phase space to 
study the dynamics. Notice that if the medium is infinite, there exists a perfect symmetry under 
translation along the array. With this in mind, we define two auxiliary variables ^ and r], that 
have sense for large enough A'^ and when the front is far enough from the boundaries: 



^ = E + Vj "^od l2^2b{r - 1)) (5) 

1 ^ 



Variable { accounts for the propagation of the front and is not bounded by default, then it 



must be defined within the range A,^ = [2yb{r — 1) + 2yb{r — l)]/y/2 = 2y26(r — 1), that is 
the magnitude that ^ increases or decreases when the front moves one position along the array. 
On the other hand, rj is defined so that it is bounded under front propagation, notice that it is 
defined from Eq. (4). If the medium consists of a number N even of units, the static solutions 
are located in this cyclic phase space in: 

A:^^ri^O (7) 



B:^ = ±j2b{r-1),7]^0 (8) 



If is odd both solutions exchange their coordinates. In what follows, we only take even, 
unless it is specified. 

The transition shown in Fig. 3 is now shown in the cylindrical phase space in Fig. 5. Notice 
that when the front is static, we have an equilibrium point at ^ = = (A-state). When 
the front oscillates a limit cycle exists, and finally when the coupling goes beyond D^h we find 
two symmetry related limit cycles, that turn around the cylindrical phase space in opposite 
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Fig. 5. Evolution on the reduced cylindrical phase space spanned by ^ and rj for three different values 
of D and r = 14; from left to right: standing, oscillating, and travelling front. At D = Dth ~ 39.63 a 
double homoclinic connection to the B-state arises. 

directions. The situation is quite similar to a pendulum where libration corresponds now to 
oscillation, and rotation of the pendulum corresponds to propagation. Thus, at D = Dth there 
exist two homoclinic loops that connect the S-state with itself (called separatrices for the 
pendulum) . This bifurcation is called gluing bifurcation because it involves the collision of two 
cycles to create another; although in our case it could be more appropriate to talk of a inverse 
gluing or a splitting bifurcation. 

Nonetheless, it is better to visuahze our gluing bifurcation in IR^. A possible transformation 
consists, topologically, in what follows: First of all, cut the cyhnder with two planes perpendic- 
ular to its axis in such a way that the saddle point and the heteroclinic orbits stand between 
both planes. Then compress both circumferences that limit the "cylinder" inbetween to a point. 
Thus, we have got an object topologically equal to a sphere. These two steps could be sub- 
stituted by making the section of the cylinder equal to zero at 77 = ±oo and then making a 
transformation that makes our object finite. The last step is to make a hole at ^ = 77 = and 
deform what remains into a plane. 

In Fig. 6 the sketch of a gluing bifurcation equivalent to the one shown in the cylindrical phase 
space (Fig. 5) is depicted. Notice that for D < Dth there exist an orbit that approaches twice 
per period close to the saddle point (the S-solution), whereas for D > Dth two symmetry 
related cycles, corresponding to both senses of propagation, coexist. 

4-2 Velocity of the front as a function of D — Dth 

The collision (and subsequent destruction) of a periodic orbit with a saddle, called saddle-loop 
bifurcation, is characterized by a logarithmic lengthening of the period of the cycle [24]. For 
the scenario depicted in Fig. 6, we expect to find the following dependencies for the periods 
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Fig. 6. Schematic of a (inverse) gluing bifurcation on M?. For D < Dth the trajectory approaches twice 
per cycle the saddle point and the front oscillates. Beyond the critical value D = Dt^ two symmetry 
related orbits, corresponding to both senses of propagation of the front, appear. and stand for 
the eigenvalues of the saddle fixed point at homoclinicity. 



(Ti 2) in a neighbourhood of Dth'- 



T, = a,-^HDth-D) (9) 
- = T2 = a2-T^ln(D-A/.), (10) 

C Au 

where Xu is the unstable eigenvalue of the saddle fixed point. Taking into account that one 
turn around the cylinder is equivalent to the movement of the front in one cell, it is clear that 
T2 is the inverse of the velocity of the front (c). In the equation for the period Ti a factor 2 
appears because the orbit approaches twice per cycle to the neighbourhood of the saddle point. 
Moreover, it is expected that the 'fast dynamics' (the motion far from the saddle) contained 
in variables ai 2 will be approximately the same at both sides of the transition, and then, 
ai ~ 2a2- In Ref. [18] it was shown that at r = 14, the behaviour of Ti/2 and agrees with 
the predictions given by Eqs. (9,10), and we show it here for r = 16 in Fig. 7. 

Notice that the velocity of the front grows quite abruptly from zero at D = Dth- The derivative 
of c = c{D), near Dth is, according to Eq. (10): 



dc if 1 



dD K{D - Dth) V«2 - ^ HD - D 



^ 00. 11, 



'thj 
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Fig. 7. 1/c (solid line) and T\j1 (squares) as a function oi D — Dth (a) and \a{\D — Dth\) (b) for r = 16. 
The behaviour agrees with that predicted by Eqs. (9) and (10). The critical coupling is found to be 
^th^^ = 20.13952. In (b) a straight line with slope — is depicted, being A„ = 1.7959 the unstable 
eigenvalue of the 5-state computed numerically for an array of A/" = 21 units. 



4-3 Quantitative analysis 



In order to verify, not only qualitatively, but quantitatively the tendency of Ti and c when 
D Dfh we computed numerically the value of the eigenvalues of the S-state at D = Dth for 
finite N. To preserve the symmetry, when computing the eigenvalues, the number of units 
was chosen to be odd, with the central unit located at the origin. The jacobian matrix (J) for 
a discrete reaction-diffusion model (see Eq. (1)) exhibits the following tridiagonal structure: 



J 



^Hi D 
D H2 D 
D H3 



V 









\ 






Hiv-i D 
D Hjvy 



;i2) 



which, according to Eq. (2) has got the components: 
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Fig. 8. Eigenvalues of the B-state for different values of D and r. 



recall that a;(jv+i)/2 = y{N+i)/2 = 0. Therefore, —6 is always an eigenvalue with trivial eigenvector 
(0, 0, 0, 1, 0, 0)"^. Using the program matlab we obtained the eigenvalues for an array of 
= 21 units, that we expect to be large enough to inform us about the properties of a front 
in an infinite medium, as long as the front is a very localized structure. 

Figure 8 shows the results for three values of the parameters r and D. Fig. 8 (a) shows the 
eigenvalues (Aj) loi D = and r = 14, which illustrates the meaning of the different eigenvalues. 
There are three simple eigenvalues corresponding to the oscillator located in the center of the 
array whose coordinates are X(Ar+i)/2 = y{N+i)/2 = Z(^n+i)/2 = 0, recall that we are considering 
a 5-state. Moreover, there are three (A^ — l)-degenerated eigenvalues, corresponding to the 
oscillators with coordinates C±. Eigenvalues are depicted with circles of different sizes to better 
observe degeneracy. When D grows above zero degeneracy is lost. In Fig. 8 (b), we graph the 
eigenvalues for D — Dt^ Ri 39.63 (only those eigenvalues with Re{Xi) > —10 are shown). The 
unstable eigenvalue corresponds to the fixed point at the origin of a Lorenz oscillator. There is 
also a slightly negative eigenvalue (As) that is the leading eigenvalue among the whole spectrum 
of negative eigenvalues. 

A different scenario is found for r = 16, as shown in Fig. 8 (c). In this case, it is not so clear 
which stable eigenvalue must be considered as the leading one in the gluing process. There are 
several eigenvalues making up two lines at both sides of the real axis in the complex plain, 
but these are not relevant, as long as they are the eigenvalues associated with the stability 
properties of each domain at both sides of the front where quasi-homogeneous domains around 
C± exist. As well, if the eigenvalue located at —b is considered, its eigenvector becomes trivial 
and it cannot be expected to participate in the gluing bifurcation. Then, the solitary eigenvalue 
denoted by \s{= p+iuj) is the one (with its complex conjugate A*) that defines the leading (two- 
dimensional) stable manifold. In fact, by increasing continuously r along Dth, the eigenvalue 
As could be followed directly from that obtained for r = 14 (see Figs. 8 (b,c)). According to 
the statements exposed above, the graph of 77 vs. ^, for D close to Dt^ and r = 16, shows some 
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Fig. 9. Projection onto (^,^7) at r = 16, D = 20.136 < Dth ~ 20.1395. In the inset, the spirahng 
approach of the trajectory to the S-state (fined circle) may be observed. 

spiralingQ when the trajectory approaches the i?-state. Figure 9 shows this effect that confirms 
that now the gluing is mediated by a saddle-focus. The cyclic definition of ^ has been relaxed 
in order to get a better observation of the spiral trajectory near the 5-state. Nonetheless, if 
one looks at the asymptotic properties of the oscillation period of the front and its velocity, 
no difference with the case of the planar connection is found. So, Fig. 7 already confirmed the 
Eqs. (9) and (10), and the slopes agree with the numerical value: ~ 1.7959. 

A sketch of the gluing bifurcation in the saddle-focus case is shown in Fig. 10. Homoclinic orbits, 
aX D = Dth are denoted by Fq and Fi. The process of gluing is drawn in analogy to Fig. 6. 
Although, from both figures, one could think that the complex and the real cases are almost 
equivalent, there exist very fundamental differences between both cases. Thus, the transition 
from oscillating to travelling fronts may become quite convoluted in the complex case. This is 
observed for larger values of r and is the subject of the next section. 



^ One should define a new extra variable C, to work in a hyper-cylindrical phase space. One could 
define, for instance, Q = '^Xji/j — bzj. Nonetheless, we shall continue to work with ^ and ij only, 
although keeping in mind that we are projecting the third variable. 
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Fig. 10. Schematic of a gluing bifurcation mediated by a saddle-focus. Two homoclinic trajectories 
Fq^i exist aX D = Dth- 

5 Exotic front dynamics 



In Figure 11 the value of the complex part of the stable leading eigenvalues (cj), along the line 



D = Dth , is shown as a function of r with a dashed line. It is found that above r 



sf 



15.45 



the stable leading eigenvalues are complex conjugates, and therefore, the gluing bifurcation 
occurs mediated by a saddle-focus point. 

Many works have been devoted to the problem of saddle-loop bifurcations involving a saddle- 
focus fixed point (see e.g. [24,27,28] and references therein). A deep presentation of this problem 
is out of the scope of this paper, so we just introduce some of the most relevant results. 



For the case of a single homoclinic orbit, occurring at a critical value of the control parameter 
(e.g. /i = 0), it was found by Shil'nikov [26] that when the ratio 



A,, 



_P_ 

A,, 



(14) 



called saddle index, is less than one, there are a countable infinity of periodic orbits in a 
neighbourhood of the homoclinic orbit, all of which are saddles. These saddle orbits are created 
in a sequence of tangent bifurcations at both sides of /i = 0, such that the stable branches 
(provided 6 > 1/2) become unstable through period-doubling bifurcations. This scenario, where 
infinite orbits, with periods ranging from some finite value to an infinite one, arrange in a wiggly 
curve around the critical value, is known as Shil'nikov wiggle. And it is, with the Lorenz and 
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Fig. 11. Complex part of the stable leading eigenvalue uj (dashed line) of the 5-state and saddle index 
5 (solid line), computed along D = Dth as a function of r. For r > Vgf, a; > 0; and for r > ri, S < 1. 



the bifocal mechanisms, one of the universal routes to chaos through homoclinicity [27,28] . 

However, the results by Shil'nikov apply in a very small neighbourhood of the critical parameter. 
Therefore, the relationship between the local (theoretical) results and the global (observed) 
behaviours cannot be stated ovcrsimply. It was shown by Glendinning and Sparrow [29] that 
the difference between the 6 < 1 and the 6 > 1 approaches to homoclinicity may not be easily 
observed numerically and may not be relevant from a global point of view. In fact, for the lowest 
branch (the orbit with the lowest period) that already exist far from the critical parameter, 
deviations from the behaviour predicted by the local analysis are most likely to occur. 

One may realize the complexity of the problem if one notices that, besides the Shil'nikov 
wiggle, there exist subsidiary homoclinic connections. Following the nomenclature of Ref. [29], 
subsidiary homoclinic orbits are multiple-pulse homoclinic orbits, i.e. orbits that pass several 
times near the stationary point without achieving homoclinicity. 

In case the homoclinic connection is double (two loops), as occurs in the gluing bifurcation, 
some differences are appreciated with the case described above. The most important point is 
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Fig. 12. 1/c (solid line) and ri/2 (squares) as a function of D fov r = 20. Fig. (a) shows that the 
oscillation period does not grow continuously. An inset, showing one of this jumps can be seen in (b). 
The velocity of the front -Fig. (d) and inset (c)- is not defined in some regions where chaotic states 
are found. Period-doublings are marked with circles and different states are labelled according to its 
symbolic code. Because of the symmetry all O's may be substituted by I's. 



that, although for 5 > 1 no Smale horseshoe exist as in the single case, the approach of orbits 
to Fo U A (Fig. 10) is chaotic [25,28]. 

Figure 11 shows the value of 6 (solid line) as a function of r along the line D = Dt^. It is found 
that at r = ri ^ 22.3, 5 crosses one. Interested in the complex behaviours that our system 
could show, we focused in two values of r above and below ri: r = 20 (5 ~ 1.18) and r = 23 
{5 0.96). We do not want to carry out a detailed analysis of the solutions found, but just to 
show that the results have got the "flavour" expected according to previous theoretical studies 
on bi-homochnicity to a saddle-focus [25,30,31]. 
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i ^ 

Fig. 13. Chaotic motion of the front for r = 20 and D = 17.85 during 200 t.u. after a transient of 
30000 t.u. a) X variable in gray scale where abcissa corresponds to the oscillator number and time is 
running downwards, b) Cyhndrical variables ^ and rj during the 200 t.u. shown in the left figure. The 
interval chosen for variable ^ has been shifted half a period, with respect to previous figures, to better 
observe the homoclinic chaos that is organized at the B-state. 

5.1 r^2Q(5>l) 

In the following, we analyse the transition from a oscillating front to a travelling front for 
r = 20. Fig. 12 (a) show the dependence of the period of oscillation of the front (Ti) vs. D. It 
is found that the period docs not grow continuously, and instead, there are some jumps, which 
give rise to hysteretic cycles (Fig. 12 (b) shows an inset of 12 (a)). These wiggles (note that 
unstable branches join consecutive stable branches) are not surprising according to [29], and 
are expected to appear w^hen approaching 5 = 1. Unfortunately, we have not been able to find 
oscillations with semi-periods larger than 3.2 t.u. This occurs because as the period grows, and 
therefore (bi-)homoclinicity is reached, the chaotic transients become larger and larger. 

The different regimes may be labelled by their representative symbolic sequences of O's and I's, 
depending on the number of turns at each side of the saddle point. For example, the uniform 
propagating regimes are labeled by {0} or {1}, but if a period-doubling occurs the new (stable) 
regimes are labeled by {0^} or {1^}, respectively. On the other hand, the oscillating regime, 
described by a two-lobed cycle, is represented by the code {01} (or {10}). 

For the travelling front region, the dependence of the velocity with D shows some features that 
were also reported in [29] , where a system of three ordinary differential equations, that exhibits 
a single homochnic connection to a saddle-focus, was studied. In Fig. 12 (c,d) the inverse of 
the velocity is plotted as a function of D. In those intervals of D where no line appear, the front 
continues to exhibit chaotic motion, characterized by "spontaneous" front reversals, after, at 
least, 30000 t.u. of transient; but some small 'windows' with regular motion can be found. 

Figure 13 shows the dynamics of the front for D = 17.85 during 200 t.u., after some transient. 
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Fig. 14. From right to left a period-doubling cascade leading to chaos. Values of D belong to the 
interval shown in the inset of Fig. 12 (c). 



In Fig. 13 (a), the variable x of the oscillators of the array is represented in gray scale, whereas 
Fig. 13 (b) shows the phase portrait rf). The front reverses its propagation several times and 
displays a quasi-erratic motion. Although the behaviour is likely to be a chaotic transient, from 
a practical point of view, it is indistinguishable of "true chaos" . 

On the other hand, transitions between different regimes exhibiting sustained propagation, 
occur at period-doubhng bifurcations (circles in Figs. 12 (c,d)). In addition, we find the well- 
known period- doubling route to chaos, see the inset of Fig. 12 (c). In Fig. 14 three periodic 
orbits and a chaotic one are shown. From right to left, a period-doubling cascade leads to a type 
of chaos that is characterized by sustained propagation of the front with non-periodic velocity. 
This scenario, that generates chaos with 5 > 1 at homochnicity, was already found in [29]. 



5.2 r = 23 f(5 < i; 



For r = 23, 5 is less than one, and therefore in a small neighbourhood of Dfu, one expects 
to find, according to [30], Shil'nikov wiggles for both, the two-lobed cycle (oscillating front) 
and the one-lobed cycle (travelling front). Also we expect to find chaos of the type reported by 
Arneodo el al. [31]; where the strange attractor organizes around a saddle-focus with symmetry, 
like in Fig. 13. 

However, we have seen in the previous subsection, that the local analysis is not always repre- 
sentative enough of the observable results that usually correspond to the lowest branch (i.e. 
the smallest period). In fact, our results do not differ too much from those found for r = 20. 
In Fig. 15, oscillating, as well as, travelling solutions are arranged according to their periods 
(considering the hyper-cylindrical phase space). 

The main oscillating solution (code {01}) disappears in a saddle-node bifurcation at ~ 18.04. 
Surely, oscillating solutions with {01} code and larger period exist, but we have not been able to 
find them. Besides the standard oscillating solution, we find others labelled {0^1^} and {0^1^} 
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Fig. 15. 1/c (solid line) and Ti/2 (squares) as a function of D for r = 23. Two periodic states with 
codes 0^1^ (dots) and 0^1^ (diamonds), as well as a conjectured unstable travelling solution (dashed 
line) are shown. 



that correspond to more sophisticated oscillating regimes. Their dynamics onto (^,77) is depicted 
in Fig. 16. Their existence can be attributed to subsidiary homochnic connections. 

Concerning the propagating solutions (solid line in Fig. 15), the "period-one" propagating 
solutions (codes {0} and {1}) undergo a period-doubling bifurcation aX D ^ 21.65. Nonetheless, 
this "period-one" solution reappears at D ~ 19.43, so we conjecture, according to previous 
studies [30], that there exist an unstable solution (dashed line) linking both regions. The lowest 
values of D with propagating solution exhibit another period-doubling bifurcation at D ~ 18.37. 

For those intervals of D, where no periodic orbits were found, chaotic behaviour (like that of 
Fig. 13) is found. Nonetheless, some very small periodic windows are found intermingled into 
the intervals with chaotic motion. 
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Fig. 16. Periodic orbits 0^1^ and 0^1*^ for r = 23 and two values of D. The cyclic restriction of the 
variable ^ has been suppressed to better observe the periodic dynamics. 



6 Further remarks 



6.1 Large D 



As the couphng D is increased, a larger amount of oscillators constitute the steepest part of the 
front. In some sense, the front becomes more 'continuous' (or less steep). This tendency allows 
us to understand better the behaviour of Dos and Dth in the large D region. When D is large, 
neighbouring oscillators have similar {x, y, z) values. Therefore, once the D^s line is crossed, a 
small increase of D (in comparison with Dth or Dos) is needed to achieve the multiple collision 
of cycles, i.e. the Dth line. As long as D is large, both A- and S-states are in a quasi-continuous 
and should exhibit similar eigenvalues spectra. 

So, as r — > Too ~ 13.5, A„ and of the i?-state at Dt^ should meet the pair of complex 
conjugates eigenvalues that characterize the Hopf bifurcation of the A-state at Dgs- In short, 
both hues meet at D = oo in a double-zero eigenvalue (Takens-Bogdanov) point. For this 
reason, oscillations above Dos have a small frequency and 5 — > 1"*" when r — > Too (see Fig. 11). 



6.2 The effect of parameter mismatch and asymmetry 



Our system of coupled Lorenz oscillators is nothing but a mathematical model, because it is 
clear that it is impossible to build an experimental set-up infinitely large with identical units 
all of them being exactly symmetric. This means that in a real system no invariance under 
translation and under reflection exists. 

Boundary effects can be significatively minimized by using an array sufficiently large, recall 
that the front is very localized, and only for large D a significative amount of oscillators are 
located in the steepest part of the front. The effect of using non-identical units is difficult to be 
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predicted but could manifest in a neighbourhood of Dth by letting some localized oscillation to 
exist or making the front to reverse back at a given time. Nonetheless, there arc lots of ways 
to make the oscillators non-identical and then giving rise to an enormous variety of effects. 

More tractable is the case where all the oscillators are identical, but asymmetric. In this situa- 
tion, we should use the theory of imperfect gluing bifurcations. Larger is the asymmetry, larger 
is the interval of D at D ~ Dth where new periodic orbits appear. Mainly, the case when the 
saddle value a = Xu + Re{Xs) is negative has been considered [33-35] showing that no more 
than two periodic orbits can coexist in the same parameter domain, and providing some simple 
rules for the symbolic sequences of O's and I's of the cycles, based on a Farey tree structure. 
Experimental studies can be found in Refs. [23] and [36]. The latter also considers the case with 
a > (i.e. saddle-focus with S < 1). 



6.3 r>rH 



The behaviour of the system when the line r — th is crossed, and C± become unstable, is 
similar to that described for the off-diagonal coupling in [32]. The transition, between oscillating 
and travelling fronts, becomes a transition between two well-differentiated types of spatio- 
temporal chaos. The slightly unstable nature of C±. makes possible that domains with oscillators 
close to C± exist, however, sometimes one oscillator is able to jump to the other lobe, which, 
in the propagating region, produces the creation of two counterpropagating fronts. In short, 
for D > Dth, we observe spatio-temporal chaos characterized by the spontaneous creation of 
counterpropagating fronts and front reversals. 



6.4 Universality? 



One may ask if the symmetry front bifurcation presented here is universal or not; in other 
words, in which kind of system could one expect to find this transition? 

It was shown in [18,32] that other coupling matrices were able to induce front propagation in 
an array of coupled Lorenz oscillators. If one considers coupling matrices with all elements zero 
except one, these two off-diagonal couplings: 



r = 7w = 4i5«2 (15) 

r = 7w = 42^/1 (16) 

exhibit front propagation by a similar route to the one explained here. However, it may happen 

that the transition occurs in such a way that the standing front loses its symmetry in a small 
interval. This happens because a pitchfork bifurcation renders the A-state unstable, and when 
the coupling in increased further, the new (static nonsymmetric) solutions "collide" with the 
S-state transferring the stability (through a pitchfork bifurcation again). Later, it is the B- 
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state which undergoes a Hopf bifurcation, and finally, the oscillating solution touch the A-state 
creating the travelling solutions. The logarithmic laws, Eqs. (9) and (10), are also obtained [18]. 

We have also checked our results with other two bistable systems: an array whose local dynamics 
is a truncation of the magnetohydro dynamic partial differential equations of a disc dynamo [37], 
and the FitzHugh-Nagumo model [9]. 

The dynamo model is similar to the Lorenz system and in some parameter range it exhibits 
symmetric bistability. The equations of the model are: 

X — a{y — x) 

y^zx-y (17) 
z — P — xy — Kz 

The system has one unstable fixed point at Pq = {0,0, (3/ k.) and two stable fixed points P± = 
(±a//3 — K, ±a//9 — K, 1) in the range P G {k, Ph), with (3h — 15 for a — 5, k — 1. We have found 
a transition like that with Lorenz systems, for instance for the coupling matrix F = 7^; = 5^2^12- 
Depending on the value of (3 the transition is smooth or chaotic. Whereas for (3 = 6 wc find a 
fine logarithmic profile of the front velocity, for (3 = 14 the front velocity function is interrupted 
by a chaotic regime, like that shown in Fig. 13, when approaching the threshold. Also the 
oscillating dynamics, for /? close to Ph is not so simple like that shown in Fig. 3. 

For the discrete FitzHugh-Nagumo model: 
iij — Uj — — Vj + D{uj+i + Uj-i — 2uj) 

Vj = e{uj - aiVj - oo) j ^1,...,N (18) 

we take Oq = and Oi = 2 which provide the Z2 symmetry and bistability, respectively. The local 
dynamics presents a sa ddle equili brium point at the origin, and two odd symmetric stable fixed 
points {u±,v±) — (±y^^i, i^y^^i). As occurs in the continuous version [16], propagation 
only succeeds for small e (0(10"''^)). For very small D, only one stable solution exits, the standing 
one (A-state). When D increases, two (counterpropagating) travelling solutions coexist with 
the standing one. Finally, the standing solution undergoes a subcritical Hopf bifurcation and 
the travelling solutions become the only stable ones. This route is apparently very different 
from those shown above, since a computation of the eigenvalues corresponding to the S-state 
reveals that A„ > — {S < 1), and therefore, if any gluing bifurcation exists, it involves unstable 
cycles. Therefore, we speculate that the transition is as follows: When D reaches a critical value 
Dsn, two stable travelling solutions (both propagation senses) are born with nonzero velocity 
{c — Co ^ 0), in two simultaneous saddle- node bifurcations. The unstable travelling solutions 
that appear in these saddle-node bifurcations become glued, when D is slightly increased, at 
D = Dg- In this way, the (unstable) oscillating solution, that coalesces with the A-state at the 
(subcritical) Hopf bifurcation, is created. 

In Fig. 17, we present the values of 1/c as a function of D, for e = 0.1. One may observe that 
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Fig. 17. Velocity of the front for the discrete FitzHugh-Nagumo model (a). Internal parameters are 
ao=0, oi = 2, and e = 0.1 (see Eq. (18)). Data arc shown with dots and circles. It may be observed that 
the travelling solution ceases to exist at D = Dsn through a saddle-node bifurcation. Nonetheless, 
in the interval shown with circles the increase of 1/c looks logarithmical. In Fig. (b) a semilog plot 
shows that, at some interval, the slope of 1/c agrees with the value of the unstable eigenvalue of 
the S-solution at Dg- A„ = 0.6781. This suggests that the travelling solution is approaching the 
(non-complex) saddle point (S-state) when D decreases. However the collision is prohibited, because 
the travelling solution is stable and the saddle index of the S-solution is approximately 0.7 for the 
values of D considered, which implies that the gluing bifurcation may occur between unstable cycles 
only. Hence, the gluing mechanism involves the unstable travelling solutions created at -D = Dsn- In 
Fig. (c) a log- log plot shows the dependence of 1/c as a function of D — Dsn- The slope near Dsn 
agrees with the expected value 0.5, characteristic of a saddle-node bifurcation. 



within a range of values of Dof D, 1/c exhibits a logarithmic profile (Fig. 17 (b)) but finally 
departs from that tendency, and shows a root-square dependence, typical of a saddle-node 
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bifurcation (Fig. 17 (c)): 



(19) 



In our case, the + (resp. -) sign corresponds to the unstable (resp. stable) solution. Dsn is 

close to Dq and that is the reason why the partial logarithmic dependence can be recognized 
in Fig. 17 (b). In fact the saddle index is not very small, 5 ~ 0.7, what makes possible the 
unstable cycles to disappear at Dsn close to Dq- 



7 Conclusions 

It has been shown that an array composed of symmetric bistable units undergoes a front 
bifurcation originating two counterpropagating travelling solutions. The mechanism consists in 
a Hopf bifurcation of the static front, followed by a global bifurcation equivalent to a gluing 
bifurcation of cycles onto a cylindrical phase space. Accordingly, close to the threshold, the 
period of oscillation of the front (Ti) and the speed of the front (c) follow logarithmic laws. An 
unstable static solution mediates the gluing process, in such a way that the value of its unstable 
eigenvalue determines the rate of divergence of Ti and c~^. The transition is typically discrete, 
and the question of the continuum limit is intriguing. 

Also, it has been demonstrated that the gluing transition may be mediated by a saddle-focus 
point. In that case, when the saddle index (5) approaches one, the transition becomes much 
more convoluted. Different oscillating and travelling regimes are observed, included chaotic 
motion of the front due to Shil'nikov chaos. 

In Sec. 6.4 we have dealt with other couplings and other bistable systems. We have found 
common features with the transition described here in detail for the Lorenz oscillator, coupled 
through the y variable. Even in the case of the discrete FitzHugh-Nagumo model, where the 
transition is apparently very different, a careful examination indicates that a similar scenario 
exists, despite the gluing bifurcation is "invisible" . 

A further investigation is needed to know a priori what local dynamics and what couplings are 
suitable to achieve propagation in a discrete symmetric bistable medium. Besides the multi- 
variable nature of the local dynamics, some trivial considerations about the coupling, and 
the observation that foci are more adequate than sinks, additional arguments are still to be 
developed. We hope that future research will clarify this point. 
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